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On the Determination of the Equations for the Perturbation of the 
Inclination and Node of an Orbit by the Method of Variation of 
Constants . By Robert Bryant, D.Sc., B.A. 


The method adopted by Oppolzer and by Watson for deter¬ 
mining the differential equations for the change in the position of 
the plane of the orbit seems to me less simple than the following. 

Taking the sun’s centre as origin, let three rectangular axes 
be chosen, whose positions will be more closely specified here¬ 
after. 

Let 

x, y, z be the coordinates of the disturbed planet, 

and 

x,y\z’ n „ disturbing „ 

Then 1c being the Gaussian constant and p the semi-parameter 
of the orbit, we have 

kVp (i +m) cos i = xv v —yv x 
7c\/p{i + m) sin S3 sin i — yv z — zv y 
— 7c Vp (1 + m) cos S3 sin i = zv x —xv z 


where v x , v y , v x , are the velocities of the disturbed planet parallel 
to the coordinate axes. 

By the differentiation of these equations and the introduction 
of suitable modifications we obtain the required results. 

Now the quantities on the right-hand sides of the above 
equations are subject to two distinct kinds of variation : one 
arising from a change in the time only, such as occurs in the 
case of undisturbed motion; the other, that arising from the per¬ 


turbation. 


d 

Denote the former of these by the symbol and the 

at 


d d 

latter bv and let denote the total variation,so that 
J dt dt 


d _d x + c? 2 
dt dt dt 


Then putting for brevity 

P -7c */p (1 + mj, 


we have 


dP . -d . .di dx , dv„ dy dv r 
— cos «-P sine . = —- v v + x —H- f-v x -y — 
dt dt dt J dt dt * dt 

^ sin S3 sin i + P cos S3 sin i + P sin S3 cos i ^ 
dt dt dt 

dy dv 0 dz dv„ 

dt " * dt dt y dt 

~ cos S3 sin i— P sin S3 sin i + P cos S3 cos i ^ 
dt dt dt 

dz , dv x dx dv £ 

~di Vx IP dt Vz x ~dt' 
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The total differentials 

dx dy dz 
dt ’ dt’ dt 

are the velocities parallel to the coordinate axes. 
Also 


at * dt \ dt dt) b \ dt dt ) 

= x d fy~ v d ^2->--- d - v ” - .• 


-y JLzZ + x. 

U dt (it 


Now at any instant the body is supposed to be describing an un¬ 
disturbed orbit, in which case the area swept out by the radius 
vector is proportional to the first power of the time only, and its 
rate of change is therefore constant, as is also the projection 
upon the coordinate planes. 

But 


X — y 

dx dt 


is proportional to the second differential with respect to the time 
of this area projected on the plane of xy , and therefore vanishes, 
as also do the corresponding quantities for the planes of yz and zx. 

Also 

^ (velocity in any direction) = acceleration in that direction due 
dt 

to perturbation. 

Hence if X, Y, Z be the perturbative accelerations parallel to 
the respective coordinate axes, our equations are simplified to 

~ cos i — P sin i~ = xY — yX 
dt dt y 

~ sin S3 sin i + P cos S3 sin i + P sin S3 cos i ~ = yZ—zY 
dt dt dt 

— cos S3 sin i — P sin S3 sin i + P cos S3 cos i~=x Z — zX. 
dt dt dt 


From the last two equations we deduce 

sin i + P cos i ^ = (yZ—zY) sin S3 + (xZ—zX') cos S3, 

CLo Cut 


whence by aid of the first equation 

P^i — (yZ— zY) sin S3 cos i + (xZ — zX ) cos S3 cos i — {xY—yX) sin *. 
dt 

Also 

P sin i ^ = (yZ—zY) cos S3 —(xZ—zX) sin S3. 
dt 

Hitherto we have left the positions of the rectangular coor¬ 
dinate axes quite arbitrary, except that their point of intersection 
is at the Sun’s centre. Our equations will be still further sim- 
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plified by a suitable choice of axes. First, let the axis of x be 
directed towards the ascending node of the instantaneous orbit 
upon the ecliptic. 

Then 


and 


P — = (xZ — zX) cos i — (xY—yX) sin i, 
dt 


P sini^ ^yZ—zY. 

at 


Next, if we adopt the plane of the instantaneous orbit as the 
plane of xy, we have still more simply 

dt 


P [ sin ‘f] =yZ - 


Sin i cZ& is the projection of the perturbation of the node 
upon a great circle perpendicular to the plrne of the instan¬ 
taneous orbit. 1 . * t ^ 

There now only remains to determine the disturbing acceler¬ 
ation Z ; X and Y now no longer appearing in our equations. 
The fundamental equations of notion referred to axes origi¬ 


nating in the Sun’s centre are 


( ~ Z \ — Z ' 

V P s / r' 3 ’ 


p being the distance between the disturbing and disturbed planets 
with two similar equations for x and y . 

We will consider only one disturbing planet, the results for 
several planets merely requiring to be added together. 

Then 

Z = k 2 m' -— —. 


But as we have adopted the instantaneous plane of the orbit 
as the plane of xy, we have 

z = o. 

Also if u be the argument of latitude, 
x — r cos u, y = r sin u, 
and substituting for P its value 


we have finally 


k */p (i + lit). 




di _ 

dt Vp (i + m) 

. . dQ> EJcm' 
sin l -= — rz=r. 


z'r cos u y 


where 


dt Vp (i + 771) 

K= i-i. 

pi r 3 


z r sm a. 
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We might also have determined from the above equations the 
value of 


and thence 


dP 
dt 9 

dp 

dt’ 


and thus have found a relation between the variations of the 
major axis and the eccentricity. Thus our equations really de¬ 
termine the variations of three of the six elements of the orbit. 


On the Variation of Latitude by the Greenwich Transit-Circle 
Observations. By S. C. Chandler. 

The interesting communication of Messrs. Thackeray and 
Turner in the Monthly Notices for November 1892 was manifestly 
prepared before seeing the final law of the variation of latitude 
in No. 277 of the Astronomical Journal. I trust they will 
not deem it an encroachment on their province to point out 
that the general accordance of the observations with the provi¬ 
sional law of No. 267, which they have used, is still more 
satisfactory when the comparison is made with the definitive 
' law. Without presuming to make a refined discussion, for which 
astronomers will look with interest to the continuation of the 
investigation by Messrs Thackeray and Turner, I find that the 
data of their Plate I. may be very satisfactorily represented by 
the expression 

+ o r,, 09 — o"*i4i cos {t —1885 Apr. 23) x o° 835-o"-i48 cos (0 — 300°) (A) 

The law on p. 100, Astronomical Journal , No. 277, gives 

— o"-i2o cos (£-1885 Apr. 9) xo°-8p5 — o"-iio cos (©-350°). 

The testimony of these observations is consequently entirely in 
harmony with the law in question, as will strikingly appear if 
the curve computed from equation (A) be entered upon 
Plate I. in the November Monthly Notices. Such a chart is 
given on Plate III.* The curve of equation (A) is represented 
by a dotted line; the observations by a zigzag continuous line, 
taking for the points of the latter the means of each successive 
three points, as given by Messrs. Thackeray and Turner, in order 
to eliminate, unobjectionably, the larger accidental errors. 

These results tend still further to confirm the doubt expressed 
on p. 71, Astronomical Journal, No. 273, as to whether tempera¬ 
ture-variation plays any significant part in the apparently 
anomalous results of this instrument, as has always been sup¬ 
posed. This matter has been the subject of a good deal of past 
investigation, but never with distinctly satisfactory conclusions ; 
and in view of the new light which this freshly discovered annual 

[* It was pointed out by Mr. Thackeray, but, unfortunately, not until after 
this plate had been printed off, that the theoretical curve is apparently one year 
in error. The portion for 1891 really belongs to 1890, and so on.— Secbetabies.] 
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